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ABSTRACT 

The Yukawa corrections of order 0(a ew m^ b ^/m^ v ) 1 0{a ew m^ b ^/m^ v ) and 
0(a ew m^ b s > /m^ v ) to the width of sbottom decay into lighter stop plus charged 
Higgs boson are calculated in the Minimal Supersymmetric Standard Model. These 
corrections depend on the masses of charged Higgs boson and lighter stop, and the 
parameters tan/J and fi. For favorable parameter values, the corrections decrease 
or increase the decay widths significantly. Especially for high values of tan/3(=30) 
the corrections exceed at least 10% for both b\ and 62 decay. But for low values 
of tan/?(=4,10) the corrections are small and the magnitudes are less than 10%. 
The numerical calculations also show that using the running bottom quark mass 
which includes the QCD effects and resums all high order tan /3-enhanced effects 
can improve much the convergence of the perturbation expansion. 



PACS number(s): 14.80. Cp; 14.80.Ly; 12.38.Bx 



1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM)||T|, |2|] is an attractive extension of 
the Standard Model (SM). In this Model every quark has two spin zero partners called 
squarks qi and one for each chirality eigenstate. These current eigenstates mix to 
form the mass eigenstates q~\ and q\. The third generation squarks are of special interest. 
This is mainly due to the reasons: large Yukawa couplings lead to strong mixing which 
induces large mass differences between the lighter mass eigenstate and the heavier one. 
This implies in general a very complex decay pattern of the heavier states. The dominate 
decay modes of the heavier squarks are the decays into quarks plus charginos/neutralinos, 
decays into lighter squarks plus vector bosons and decays into lighter squarks plus Higgs 
bosons. All these squark decays have been extensively discussed at the tree- level |3|, f|, M. 
The next generation of colliders, for example, CERN Large Hadron Collider (LHC) will 
be able to produce such kind of particles with masses up to 2.5 TeV0, and a e + e~ lin- 
ear Collider (LC)@ will be able to make precise measurements of their properties. Thus 
a more accurate calculation of the decay mechanisms beyond the tree-level is necessary 
to provide a solid basis for experimental analysis of observing these decays at the next 
generation of colliders. Up to now one-loop QCD and supersymmetric (SUSY) QCD cor- 
rections to the squark decays have been calculated^, and the Yukawa corrections to 
the squark decays into quarks plus charginos/neutralinos also were given in Ref.f9j. Very 
recently, a complete one-loop computation of the electroweak radiative corrections to the 
above processes has been presented by J. Guasch, W. Hollik and J. Sola[[L(J. But the elec- 
troweak radiative corrections to the heavier squark decays into lighter squarks plus vector 
bosons and decays into lighter squarks plus Higgs bosons have not been calculated yet, 
even Yukawa corrections to these processes. In this paper, we present the calculations of 
the Yukawa corrections of order 0(a ew m^/my i r), 0(a ew m^/m'y V ) and 0{a ew m^, b \/rriy V ) 
to the width of sbottom decay into lighter stop plus charged Higgs boson, i.e. the decay 
bi — > ti + H~ , where t\ is the lighter stop. These corrections are mainly induced by Yukawa 
couplings from Higgs- quark- quark couplings, Higgs-squark-squark couplings, Higgs-Higgs- 
squark-squark couplings, chargino(neutralino)-quark-squark couplings, and squark- squark- 
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squark-squark couplings. 

Our results can be generalized straightforwardly to the decay t 2 — > h + (h°, A ). As for 
the heavier squark decays into lighter squarks plus vector bosons, the electroweak radiative 
corrections are simple because of the relatively less renormalization parameters involved. 

2 Notation and tree-level result 

In order to make this paper self-contained, we first summarize our notation and present 
the relevant interaction Lagrangians of the MSSM and the tree-level decay rates for hi — > 

h + H-. 

The current eigenstates q L and q R are related to the mass eigenstates q~i and q 2 by: 

h ), R"=( ™\ sin M (i) 

with < 6q < Tv by convention. Correspondingly, the mass eigenvalues and rriq 2 (with 
m q L < m q 2 ) are given by 

with 

m\ L = M^ + m 2 q + m 2 z cos 2/3(I$ L - e q sin 2 6W), (3) 

m L = M {V,d} + m 2 q + m 2 z cos 2/3e ? sin 2 6> w , (4) 

a q = A q — yu{cot/5,tan/3} (5) 

for {up, down} type squarks. Here M| is the squark mass matrix. Mq ^ ^ and are 
soft supersymmetric breaking parameters and ji is the Higgs mixing parameter in the 
superpotential. 1\ L and e q are the third component of the weak isospin and the electric 
charge of the quark q, respectively. 

Defining H k = (h°, H°, A , G°, H ± , G ± ) (k=l...Q), one can write the relevant lagrangian 
density in the (51,52) basis as following form (i, j=l,2; a and (3 flavor indices) 

Relevant = H k f(a%P L + b%P R )q a + (G* )^5?* $ + gq{4 k P R + bi)x° k qi 
+gt{l\ k P R + k\ k P L )xt~bi + gWikPR + kl k P L )xt% 

+ (Gl k ) t] H + H k qP*q? + h.c, (6) 

2 



with 

(Gth = [R a Gt(Rh T h, (Gt k h = [R 6l G^ k (R^) T }ij (k = 1...6), (7) 

where G k and G§ k are the couplings in the (qi, Qr) basis, and their explicit forms are shown 
in Appendix A. The notations a k , b k (k=1...6), and aj k , b q ik (k=1...4), and l 9 k , k\ k (k=l,2) 
used in Eq.@ are defined also in Appendix A. 

The tree-level amplitude of hi — > t\H~ , as shown in Fig. 1(a), is given by 

M (o) _ jg ipb ( mltan(3 + cot f3 -Tripsin 2(3 m t (A t cot (3 + /i) \ ( i . T] 
* VW { m b (A b t a nf3 + fi) 2m t m 6 /sin2/5 J ^ Jil ' w 



and the decay width is 



M( 0) \ 2 \V 2 (ml,ml,m 2 H - 
167rm? 



(o) i— t i vww-h-j . v 

i 1c .3 v ^ 



6i 

with A(x, y, z) = (x — y — z) 2 — 4yz. 

3 Yukawa corrections 

The Feynman diagrams contributing to the Yukawa corrections to foj — > £i-ff~ are shown in 
Figs.l(b)-(f) and Fig.2. We carried out the calculation in the t'Hooft-Feynman gauge and 
used the dimensional reduction, which preserves supersymmetry, for regularization of the 
ultraviolet divergences in the virtual loop corrections using the on-mass-shell renormaliza- 



tion scheme [0], in which the fine-structure constant a ew and physical masses are chosen 
to be the renormalized parameters, and finite parts of the counterterms are fixed by the 
renormalization conditions. The coupling constant g is related to the input parameters e, 
raw and mz via g 2 = e 2 j s 2 w and s 2 w = 1 — m|r/m|. As for the renormalization of the 
parameters in the Higgs sector and the squark sector, it will be described below in detail. 
The relevant renormalization constants are defined as 

m wo = m w + Smly, m| = ra\ + 5m 2 z , 
m q o = m q + 5m qj m| = m| + 5m\, 
A q o = A q + 5A q , fi = fj, + 5fi, 
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0qo = 9q + 59q, tan /3 = (1 + 8Zp) tan (3, 

q i0 = (l + 5Z!) 1 / 2 + 5Zf j q j , 

H = (1 + 5Z H -) l ' 2 H- + 5Z HG G-, 

G = (1 + c^-) 1 / 2 ^- + 5 ZquH (10) 



with q = t,b. Here we introduce the mixing of and G 12], instead of the mixing of 



F- and W~ 13 



W~ = (1 + 5 Z\y-)^^W~ + iZ WH d,H-, 

H = (l + dZn-y^H-, (11) 

due to the reason: the former can successfully cancel the divergences for the case we are 
considering, and the latter, however, is just right for the renormalization of the parameters 
f3 and a, and in the case that the particles interacting with the charged Higgs boson are 
the on-shell fermions (where a /p arising from d^H~ and the vertex H^~-fermion-fermion 
is inserted between two on-shell fermions, and turned into the fermion masses by Dirac 
equations, which is therefore similar to the structure of the Yukawa coupling if~-fermion- 
fermion) . 

Taking into account the Yukawa corrections, the renormalized amplitude for bi —>■ tiH~ 
is given by 

M ren = ^(0) + gM (v) + ^ 

where 8M^ v) and 8M\ c ) are the vertex corrections and the counterterm, respectively. 
The calculations of the vertex corrections from Fig.l(b)-l(f) result in 

i 6 - 

k=l j 



i 4 - 
fc=i j,i 
ig 2 4 

— ^2{[ m tm b a b 5 + (m? - p~ b . ■ p H ~)bl]m^b\ k a%C 
° n k=i 
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+b\ k [m b a\b%(jp bi -Ph-) + m t b\b% Vli - m 2 ob b 5 a\* k (2p~ bi - p H -)U-p£C 



+f H -C l2 ) - (m b a b 5 b[* k + m t b\b% - m^b b 5 a\* k )b b k (mlC 21 + m 2 H -C 22 

-2Pfe. ■ Ph-C 23 + g£C 24 ) + (a <-> b)}(-p- bi ,p H -,m^o k ,m b ,m t ) 
3i 
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16vr 2 , 

3,1 



^(hlR^R'jiRli + h^R^R^R^G^MpH^m^m^), (13) 



where i?o and Cj(j) are two- and three-point Feynman integrals|b|], respectively, and h t ( b ) 
is the Yukawa coupling defined by 

^ = ^ = a' (14) 

\f 2mw sin p y zmjy cos P 

In the first line of 5M^ V \ q = t for k = 1...4 and q = b for k = 5, 6, respectively, while in 
the second line, q = b for = 1...4 and q — t for k = 5,6, respectively. 
The counterterm can be expressed as 

5M\ C) = i(59- b + ^(Gjbs-i,! + i(<S0 t - + *^fi)(G?i)i2 + i[R l (5Gl)(Rt) T }u 

+ 1 -{5Z\ + + 5Z H -)(Gl) a + i8Z GH {Gl) n (15) 



with 



(5G\)n = —/= [( — 5— tan f3 + vr? t cot 0) + 2m b tan /3<5m{, 

y2m w # 2 m,^ 



+2m t cot /35m t + m b 5 tan /3 + m 4 5 cot (16) 

(6G\) 12 = - + —)m t (A t cotp + fi) + m t (8A t cot (3 

V2m w g 2 77% m t 

+A t <Jcot/? + <fyx)], (17) 

(5G|) 21 = -7^— [(- - + —)m b (A b tan/5 + //) + m 6 (<L4 6 tan/3 

+A 6 6tan/3 + 5/j)], (18) 
/rA&x 2gm b m t 5g 1 Sm^ 5m b 8m t 

{oGVj22 — —?= ir~ H 1 coszpoZfl). 19 

5 y^m w sin2/3 V 5 2 m 6 m t y p ' 

Here we consider only the corrections to the Yukawa couplings. The explicit expressions of 

some renormalization constants calculated from the self-energy diagrams in Fig. 2 are given 

in Appendix B, and other renormalization constants are fixed as follows. 

For 5Zghi using the approach discussed in the two-Higgs doublet model (2HDM) in 



12 1 , we derived below its expression in the MSSM, where the version of the Higgs potential 



is different from one of Ref. |12[]. First, the one- loop renormalized two-point function is 
given by 



iF GH {p 2 ) = i(p 2 ~ m 2 H ~)5Z HG + ip 2 SZ GH - iT GH + iT, GH (p 2 
where T G h is the tadpole function, which is given by 



T, 
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GH 



■[Th 2 sin(a — (3) + T Hl cos(a — (3)]. 



2m w 

Next from the on-shell renormalization condition, we obtained 

1 



6Z, 



GH 



ni 



■[T GH - J: GH (m H - 



(20) 



(21) 



(22) 



The explicit expressions of and the tadpole counterterms T Hk {k = 1,2) are given in 
Appendix B. 

For the renormalization of the parameter (3, following the analysis of Ref. |TB[ , we fixed 
the renormalization constant by the requirement that the on-mass-shell H + fv T coupling 



remain the same form as in Eq. (3) of Ref. []T3| to all orders of perturbation theory. However, 
with introducing the mixing of H~ and G~ instead of H~ and W~ , the expression of 5Zp 
is then changed to 



5Za 



1 6m2w 1 5ml + 1 5 ™ 2 *- 5 ™} _ i 6Zh+ + cot mGH . 



2 m 



w 



2 mi 



2 mi 



m 



(23) 



For the counterterm of squark mixing angle 8q, using the same renormalized scheme as 



Ref. [Bfl , one has 



^ = Re[S? 2 K) + S? 2 (m|)] 



(24) 



2(m| - m 2 hJ 

For the renormalization of soft SUSY-breaking parameter A q , we fixed its counterterm 



by keeping the tree-level relation of A q , m qi and 6 q [15], and get the expression as following: 

+ 



,U„ = m \ m '" 2 (2 cos 26 q 56 q - sin2%^) + ^^(5m? - 6m 2 



2m n 



2m n 
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+{cot P, tan (3}5jj, + 5{cot (3, tan j3}fi. 



(25) 



As for the parameter [i, there are several schemes[10|, |16|, Ol to fix its counterterm, and 



here we use the on-shell renormalization scheme in Ref. [17], which gives 



Sfi = ^2[m^+(SU k2 V k2 + U k2 5V k2 ) + 5m^+U k2 V k2 ] 



(26) 



k=l 



where (U, V) are the two 2x2 matrices diagonalizing the chargino mass matrix, and their 
counterterms (SU, 8V) are given by 

8U = h,8Z R -5Z T R )U, (27) 
5V = ^(5Z L -5Zi)V. (28) 

The mass shifts 8m^+ and the off-diagonal wave function renormalization constants 5Z R ^ 
can be written as 

5m~+ = iRe[m~+(n£ fe (m|+) + II*(m|+) ) + n2f(mj+) + n&V#)l. (29) 

X^ Xj 

+ U ij R ( m l+) m x+ + U ij L ( m l+) m x+]i ( 30 ) 
(^x)y = (^«)ii (31) 

The explicit expressions of the chargino self-energy matrices and H S > L ( R ) are given in 

Appendix B. 

Finally, the renormalized decay width is then given by 

r< = rf } + sr\ v) + 8r\ c) (32) 



with 



X^'^(ui~ 771- TT? — ) 

5r S a) = olf H Re{Mi 0) ^Mf a) } (a = v, c). (33) 



iam- 
bi 



4 Numerical results and conclusion 

In the following we present some numerical results for the Yukawa corrections to sbottom 
decay into lighter stop plus charged Higgs boson. In our numerical calculations the SM 



parameters were taken to be a ew (mz) = 1/128.8, mw = 80.375GeV, m z = 91.1867GeV|T% 
and m t = 175.6GeV. In order to improve the convergence of the perturbation expansion, 
especially for large tan f3, we take the running mass m b (Q) evaluated by the next-to-leading 



order formula 19 



m b (Q) = U 6 (Q,m t )U 5 (mt,m b )m b (m b ), (34) 



where we have assumed that there are no other colored particles with masses between scale 
Q and m t , and mb{mb) = 4.25GeV|20[|. The evolution factor Uf is 

u,{Q„ Qi) = (^#|r [i + <^)-<^) jif\ 

sn = _J2 = .8982- 504/ + 40/^ 

33-2/' 3(33 -2/) 2 ' v ; 



where a s (Q) is given by the solutions of the two-loop renormalization group equations ||21||. 
When Q = 400GeV, the running mass rrib{Q) ~ 2.5GeV. In addition, we also improved the 



perturbation calculations by the following replacement in the tree-level couplings [19] 



where 



mm - 1 + A.MAW) <36) 



2a s h 2 
Am b = tan /5J(m^,m^ 2 , M 3 ) + j^Mt tan/W^m^, m h , /x) 

fiM 2 tan /9[cos 2 %J(mj , , M 2 , /x) + sin 2 9fl(mf 2 , M 2 , /x) 



16tt 



- cos 2 e~ b I(m~ h ,M 2 , /x) + ~ sin 2 9 b I(m~ b2 ,M 2 , /x)] (37) 



with 



J(fl ' 6 ' C) = (^-^(^-c 2 )(a 2 -c 2 ) (a2&2 bg ^ + &V bg ? + C2 ° 2 bg ^ (38) 



The two- loop leading-log relations [22] of the neutral Higgs boson masses and mixing 
angles in the MSSM were used. For m H + the tree-level formula was used. Other MSSM 
parameters were determined as follows: 

(i) For the parameters M\, M 2 and /x in the chargino and neutralino matrix, we take 
M 2 and /x as the input parameters, and then use the relation M\ = (5/3)(g' 2 /g 2 )M 2 ~ 
0.5M2@, |23| to determine M\. The gluino mass Mg in Eq.(p7|) was related to M 2 by 
M g = (a s (Mg)/a 2 )M 2 §. 

(ii) For the parameters m 2 ^ ^ ^ and A tj b in squark mass matrices, we assumed Mq = 
Mjy = Mf) and A t = A b to simplify the calculations. 

Some typical numerical results of the tree-level decay widths and the Yukawa corrections 
are given in Fig.3-9. 
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Fig. 3 and Fig.4 show the dependence of the results of b\ and b 2 decays, respectively. 
Here we take m A o = 150GeV, fi = M 2 = 400GeV and A t = Ab = ITeV. The leading terms 
of the tree level amplitude are given by 

Mf) „ [mt(A t cot (3 + n)RiR\ 2 + m 6 (A 6 tan/5 + /^i^ii], ( 39 ) 

y2m w 

where cos0 s ~ (0.54, 0.67, 0.70) and cos# t ~ ~ (-0.71, -0.71, -0.71) for tan/3 = (4, 10, 
30), respectively. In the case of % = 1, two terms in Eq.(|39"D have opposite signs, and 
their magnitudes get close with the increasing tan/3 and thus cancel to large extent for 
large tan/5. Therefore, the tree level decay widths have the feature of r (tan/5 = 4) > 
r (tan/5 = 10) > r (tan/5 = 30) in the most of parameter range, as shown in Fig.3(a). In 
the case of i — 2, for tan /3 — 4, 10 and 30 two terms in Eq. (P9"D have the same signs, so T 
is larger than the one of the case of i — 1. From Fig.3(b) and Fig.4(b) one can see that 
the relative corrections are sensitive to the value of tan /5. For tan f3 = 30, the magnitudes 
of the corrections can exceed 10% when > 160GeV for b\ decay and > 260GeV 
for b 2 decay. Especially in the case of hi decay it even can reach 40%, which is due to the 
fact that the corresponding tree-level decay width already becomes very small. There are 
the dips at =311GeV and 390GeV on the solid line in Fig.4(b), which come from the 
singularities at the threshold points m~ b2 = + m t and m~ bi = m-+ + m t , respectively. 
However, for tan/5 = 4 and 10, the corrections of two sbottom decays are always small and 
range from —5% to 5%. In general, for low tan/5 the top quark contribution is enhanced 
while for high tan f3 the bottom quark contribution become large, and for medium tan j3, 
there are not any enhanced effects from Yukawa couplings. So the corrections for tan (3 = 
4 are generally larger than those for tan/5 =10, as shown in Fig. 3(b) and Fig. 4(b). 

Fig. 5 (Fig. 6) gives the tree-level decay width and the Yukawa corrections as the functions 
of m H - in the case of 61 decay (6 2 decay). We assumed = 170GeV, \i = M 2 = 400GeV 
and A t = A b = ITeV. The features of the tree level decay widths in Fig.5(a) and Fig. 6(a) 
are similar to Fig.3(a) and Fig. 4(a), respectively. From Fig. 5(b) and Fig. 6(b) we can see 
that the relative corrections decrease or increase the decay widths depending on tan/5. 
Fig. 5(b) shows that the corrections for tan/5 = 4 are always positive and range between 
6% and 3%. For tan/5 = 10 the corrections are negligiblely small. For high tan/5(= 30) the 
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corrections exceed 10% when m#- < 180GeV. There is a dip at m^j- ~ 178GeV on the 
solid curve due to the singularity of the charged Higgs boson wave function renormalization 
constant at the threshold point m#- = m t + m b . Fig. 6(b) shows that the corrections are a 
few percent for tan/3 = 4, 10 and 30. There are a dip and a peak on the solid curve, which 
arise from the singularities at the threshold points mg = mg + m A o and m#- = m t + m&, 
respectively. 

In Fig. 7 and Fig. 8 we present the tree level decay widths and the Yukawa corrections 
as the functions of /i in the case of b± — > t% + H~ and 62 — ^ ^1 + H , respectively, assuming 
m k = 170GeV, M 2 = 400GeV, A t = A b = ITeV and m A o = 150GeV. When \i takes 
some values, the tree level decay widths are getting very small (< 10~ 3 GeV), as shown in 
Fig. 7(a) and Fig.8(a), and the corrections near the above values do not have a physical 
meaning. So we cut off the corrections, since perturbation theory breaks down there. From 
Fig. 7(a) (Fig. 8(a)) we can see that there is a high peak (a deep dip) for tan/3 = 30 and 
/i ~ 30GeV. This is due to the fact that when tan (3 — 30 and fi ~ 30GeV the second term 
in Eq. (|39|) is enhanced (suppressed) greatly for sin Q~ h ~ 1 (cos^ ~ 0) as a result of the 
off-diagonal elements of M| in Eq.(^) approaching to zero. Fig.7(b) and Fig.8(b) show that 
the Yukawa corrections depend on /1 strongly. Especially, in the region of Tq getting very 
small, the corrections get the rapid variations between the positive and negative values 
with the changes of the sign of the tree level amplitude. In general, when the tree-level 
decay widths for tan f3 = 4 and 10 are not close to zero, the corrections are always small. 
Comparing Fig. 7(b) with Fig.8(b), we can find that the Yukawa corrections for tan/5 = 30 
are more significant for b\ decay than for 62 decay. 

Finally, in Fig. 9 we show the Yukawa corrections as a function of tan/3 in three ways of 
perturbative expansion: (i) the strict on-shell scheme (the dashed line), where the bottom 
quark pole mass 4.7GeV was used, (ii) the QCD-improved scheme (the dotted line), in 
which only QCD running bottom quark mass mb{Q) in Eq.([34]) was used, and (iii) the 
improved scheme (the solid line), in which the replacement Eq.(^) was used. Here we 
assumed m h = 170GeV, fi = M 2 = 400GeV, A t = A b = ITeV and m A o = 150GeV. One 
can see that the three ways all give small corrections (|5r/r | < 5%) for tan/3 (< 15). 
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However, the magnitude of the corrections in the case (i) increases rapidly for tan/3 > 15, 
and when tan/3 > 33 the corrections will result in the physically meaningless negative 
width. But the convergences in the cases (ii) and (hi) are much better, and especially in 
the case (iii) the magnitude of the corrections still is less than 40% for high values of tan j3 
(=40). 

In conclusion, we have calculated the Yukawa corrections to the width of sbottom 
decay into lighter stop plus charged Higgs boson in MSSM. These corrections depend 
on the masses of charged Higgs boson and lighter stop, and the parameters tan/3 and 
fi. For favorable parameter values, the corrections decrease or increase the decay widths 
significantly. Especially for high values of tan/3(=30) the corrections exceed at least 10% 
for both b\ and 62 decay. But for low values of tan/3(=4,10) the corrections are small and 
the magnitudes are less than 10%. The numerical calculations also show that using the 
running bottom quark mass which includes the QCD effects and resums all high order tan /3- 
enhanced effects, as given in Ref . |H| , can improve much the convergence of the perturbation 
expansion. 
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Appendix A 



The following couplings are given in order 0(h t , hi,). 
1. squark - squark - Higgs boson 
(a) squark - squark - h° 
( 

G\ = 



V 



—\/2m q h q 



s a 



ysh q (A q 



V2 



f c a 




s a ) 








I - S <*J 





■\f2m q h 



(40) 
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for 



up 
down 



type squarks, respectively. We use the abbreviations s a = sin 
a is the mixing angle in the CP even neutral Higgs boson sector, 
(b) squark - squark - H° 

( 



G% 



— \[2m q h q 

-^h q (A q l a j 



So 
Cn 



s a 



—\/2m q h g 



(c) squark - squark - A 



( 



2m w 



V 



A n 





cot (3 
tan (3 



-A, 



cot (3 
tan/5 





- n 



(d) squark - squark - G° 



( 



Gl = 



■ grriq 

i — 

2mw 







-A q + /i 



V 



A q ~H 



cot (3 
tan/5 



cot/5 1 \ 
tan/3 J 

I 



(e) squark - squark - H^ 1 



ml tan (3 + mf cot (3 m t (A t cot (3 + //) 



y/2m w V m b (A b tan/5 + //) 2m t m fe /sin2/5 
(f) squark - squark - G ± 



t\T 
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ml — ml m t (A t — /i cot (3) 



~ V2m w V rn b (fitan(3 - A b ) 
2. quark - quark - Higgs boson 



4 = 



/if, sin /3 1 



--U 



cos (3 



y/2 q \ sin/5 J '2m 



14' 



-m 6 



^cos/5 j ' ^m^ \ m t 








' 2m 



H' 
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3. quark - squark - neutralino 



q _ p9y I 
lx i2 1 q 



ik 



iA _ _ d9 y \ N > 

°ik — rL il I q \ jy 



ki 



(4* 



v N k3 j ' I JV fc3 y 

Here N is the 4x4 unitary matrix diagonalizing the neutral gaugino-higgsino mass matrix 
0, [23] , and the Yukawa factor Y q = h q /g. 
4. quark - squark - chargino 

( V k2 1 



79 _ p9 y 
t ifc — rL i2 I q 



u. 



(49) 



y fc2 J [ -ftVfc2 J 

Here U and V are the 2x2 unitary matrices diagonalizing the charged gaugino-higgsino 
mass matrix [|2], |23| . 

5. squark - squark - Higgs boson - Higgs boson 

(a) squark - squark - H~ - (k=1...4) 

r b = (c l ) T = g2 ( m t Sk + m2bTk ° ^ 

5k { 5k) 2y/2m 2 w \ 2m t m b / sin 2(3V k J 

with 



(50) 



Sk = (cos a cos /?/ sin 2 /?, sin a cos (3/ sin 2 f3, — zcot 2 /3, zcot/3) (51) 
Tk = (— sin a sin/3/ cos 2 [3, cos a sin /?/ cos 2 /5, itan 2 /3, ztan/3) (52) 
Vfc = (sin(/3 — a), cos(/3 — a), 0, i) (53) 



(b) squark - squark — H - H + 



G q 



55 



V 



hi sin 2 /3 
/i 2 cos 2 (3 







(c) squark - squark H - G + 

( f 



56 



— m 2 tan /? 
m? cot /? 





cos 2 (3 
sin 2 /? 





cot (3 
- tan (3 



(54) 



(55) 
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Appendix B 



We define q — t and b, q' the SU(2) L partner of q, and q" = q for k = 1...4 and q" = q 
k — 5, 6. Then we have 

^ = 1fi T 2 K + m t~ M™1) - M™t) ~ m 2 t B - (m 2 - m 2 )^] 

m W 107T Tflyj 

{m 2 w ,m b ,m t ), 

5j 4- = F^V E ihlh - e q sin 2 9 W ) 2 + e 2 q sin 4 ^][2mJ - 2A (m 2 q ) - m 2 q B 
m z 57T m w q=t b 6 



-2m 2 e q sin 2 W (I% L - e q sin 2 6 w )B }(m 2 z , m q , m, 



1)1 

5Z H - = JL{2[(<4&| + fo*4)(m 2 f+ G' 1 + B! + m 2 G ) + m t m 6 (a*^ + &«&g)G ] 
(m 2 H+ , m b , m t ) - Yl G \ ) ji (^s) y G o (m 2 H+ , m~ h , m~ t . ) } , 

^ = E {2K + 6|)mAK) - E(Gf)«4.K)}> 

3 

+b\b b 6 )B ](m 2 H+ ,m b ,m t ) - J2( G ih( G e)i3 B o(™H+: m- bi , m~ t .) 

3,1 

+ E E(4)^oK)l, 

(/tq 107T j, =1 //tg Z 

4 777,-0 1 

VEEl^i^o - ^(I4 fc | 2 + |^| 2 )Bi]K,m,o,m g .) 

V E EtdH^o - ^ (l^l 2 + l4l 2 )^K, m sJ )}, 

fc=l j m 9 2 ^ 

k=l j k=l 

x(m|_Bi + A (m|o) + m 2 B ) + 2m 9 m^oRe(a^fef^)73o](m|, m„ m^o) 
-2/ EKI4I 2 + l4| 2 )K-'5i + AoK + ) + m 2 So) 



fc=l 

+2m^m~+Re(Z? fe /4*)£ ](m| , m,/, 



^ = ttL{- E E(G2)«(^')^oK,mff,m^) + 2/ £[(|4| 2 + l&? 



Kir; 2 , , 

k=l j fc=l 



LP) 
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x(£?i + m'-.Gi + m 2 Go) + 2m q m^oRe(a q k bf k )Go](m^.,m q ,m^o) 
+*9 2 EKI4I 2 + I4l 2 )(£i + m|G?! + rnjGo) 

k=l 

+2m q/ m^+Re(lf k k g i *)G }(ml,m gl ,m-+)}, 

x(p 2 B l + A)(w|o) + m 2 q B ) + m q m^a\ k b q 2k + af k b\ k )B ](p 2 , m q , m^o) 



-2 9 2 + K k kf k ){p 2 B l + A)(m 2 + ) + m JS 

fe=i fe 

+"V m v+ ( fc 2fc + ^2fc fc ?fc ) A)] (P 2 , m,' , m~+ ) }, 



Xi v 1K zlt are l«/ "J ' y ' Xfe 
1 



3 2 - - 
fc=i 

3 2 . . 



16vr 2 
3 



fc=i 

2 



n^' (p 2 ) = — -J E[ m ^fci4j 5 o(p 2 , m b, m k ) + m t l b ki k b kj B {p 2 , m t , m lk )}, 

k=l 

3 2 

RifiP 2 ) = E K44jA)(p 2 , m b, m ik ) + mk b ki l b kj B (p 2 , m t , m~ h )). 



k=l 

Here A and i?i are one- and two-point Feynman integrals [jl4j], respectively, and Gi = 
dBi/dp 2 . 
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(a) 



(6) 



H~ 

(c) 



Ax 
% 



b. 



(d) 



(/) 



Figure 1: Feynman diagrams contributing to supersymmetric electroweak corrections to 
bi — > t\H~: (a) is tree level diagram; (6) — (/) are one- loop vertex corrections. In diagram 
(6) the subscript k of H k can take from 1 to 4. In diagram (d) q — t for k — 1...4 and g = 6 
for fc = 5, 6. In diagram (e) q — b for fc = 1...4 and g = t for fc = 5, 6. 
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Figure 2: Feynman diagrams contributing to renormalization constants. In diagram 
q = t{b) for k = 1...4 and q = b(t) for k = 5,6. 
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Figure 3: The tree-level decay width (figure (a)) of b\ — > t±H and its Yukawa corrections 
(figure (&)) as functions of for tan/5 = 4, 10 and 30, respectively, assuming m^o = 
150GeV, n = M 2 = 400GeV, A t = A b = ITeV and M Q = M = M 6 . 
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Figure 4: The tree-level decay width (figure (a)) of 62 — > t±H and its Yukawa corrections 
(figure (&)) as functions of for tan/5 = 4, 10 and 30, respectively, assuming m^o = 
150GeV, n = M 2 = 400GeV, A t = A b = ITeV and M Q = M = M 6 . 
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Figure 5: The tree-level decay width (figure (a)) of b\ — > t\H and its Yukawa corrections 
(figure (b)) as functions of m H - for tan/5 = 4, 10 and 30, respectively, assuming = 
170GeV, n = M 2 = 400GeV, A t = A b = ITeV and Mq = M = M D . 
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Figure 6: The tree-level decay width (figure (a)) of 62 — > t±H and its Yukawa corrections 
(figure (b)) as functions of m H - for tan/5 = 4, 10 and 30, respectively, assuming = 
170GeV, n = M 2 = 400GeV, A t = A b = ITeV and Mq = M = M D . 
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Figure 7: The tree-level decay width (figure (a)) of b\ — > t±H and its Yukawa corrections 
(figure (b)) as functions of \i for tan (3 — 4, 10 and 30, respectively, assuming = 170GeV, 
M 2 = 400GeV, A t = A b = ITeV, m A o = 150GeV and \J () U, M b . 
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Figure 8: The tree-level decay width (figure (a)) of 62 — > t±H and its Yukawa corrections 
(figure (b)) as functions of \i for tan (3 — 4, 10 and 30, respectively, assuming = 170GeV, 
M 2 = 400GeV, A t = A b = ITeV, m A o = 150GeV and \J () U, M b . 
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Figure 9: The Yukawa corrections of 62 - ► t\H as a function of tan/?, assuming = 
170GeV, fi = M 2 = 400GeV, A t = A b = ITeV, m A o = 150GeV and Mq = M = M b . 
The dashed line corresponds to the corrections using the on-shell bottom quark mass, the 
dotted line to the improved result only using the QCD running mass mb(Q), and the solid 
line to the improved result using the replacement Eq.flSBT). 
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